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SECTION  I 


INTRODUCTION 

The  understanding  of  the  failure  mechanism  of  a  structural  material 
is  basic  to  design.  Static  design  techniques  are  assumed  to  be  adequate 
for  structural  application;  which  includes  advanced  fiber  composites.  The 
failure  surface  concept  assumes  failure  occurs  when  the  load  which 
produces  combined  stresses  falls  outside  the  defined  envelope.  The  stress 
envelope  is  produced  from  simple  static  test  results.  The  prediction  of 
fatigue  failure  using  a  static  stress  envelope  is  questionable.  Actual 
failure  is  established  by  test  of  structural  components. 

The  typical  structural  failure  will  be  induced  from  stress  fields 
created  by  discontinuities.  Most  structural  components  are  subjected  to 
cyclic  loading  during  service.  Fatigue  cracks,  originating  from  flaws 
within  the  composite,  are  propagated  to  fracture  under  cyclic  loading 
conditions.  The  fatigue  or  cyclic  loading  best  fits  structural  design. 

Three  approaches  used  in  fatigue  analysis  include:  (1)  cumulative 
damage,  (2)  finite  time  to  crack  nucleation,  and  (3)  presence  of  internal 
flaws  (Fracture  Mechanics).  Fracture  mechanics  assumes  all  materials  contain 
internal  flaws  with  unspecified  geometry  and  distribution.  The  geometry 
and  distribution  of  these  internal  flaws  must  be  established  by  NDT.  Flaws 
grow  in  length  when  the  lammate  is  subjected  to  load,  static  or  cyclic, 
until  instability  occurs.  At  this  point,  fast  crack  propagation  occurs  and 
the  end  result  is  fracture. 
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Since  the  geometry  and  distribution  of  internal  flaws  may  not  be 
determined  accurately,  a  flaw  of  known  geometry,  is  placed  in  the  test 
specimen  to  establish  failure  criteria.  Cracks  are  used  as  the  flaw 
and  the  critical  stress  intensity  factor  (K^)  determined  by  experiment. 

The  stress  concentration  at  the  tip  of  the  cracks  is  not  defined  and 
makes  reproducibility  difficult.  Slits  with  known  stress  concentrations 
at  the  tips  are  used  as  standards.  A  simple  flaw  geometry  would  be  a 
circular  hole  in  the  center  of  the  test  specimen.  The  circular  hole 
allows  calculation  of  the  stress  concentration  distribution  and  establishes 
a  basis  for  a  generalized  failure  theory.  The  stress  concentration 
distribution  around  the  periphery  of  the  hole  can  be  established  for 
static  and  cyclic  loading  conditions  and  varies  with  composite  system, 
fiber  orientation,  specimen  geometry,  diameter  of  hole,  etc.  The  maximum 
stress  concentration  should  be  the  point  of  crack  initiation  when  considering 
isotropic  materials.  The  crack  will  propagate  in  the  field  of  the  smallest 
strain  energy  intensity  gradient.  This  would  suggest  that  the  crack  can  propa¬ 
gate  in  more  than  one  direction  (also  form  more  than  one  crack) .  Crack 
instability,  however,  does  not  necessarily  mean  structural  failure. 

A  generalized  failure  theory  should  consider  as  many  design  variables 
as  possible.  These  variables  would  include  composite  system,  fiber 
orientation  (laminate  mechanical  properties),  type  of  loading,  loading 
conditions  (static  or  cyclic),  strain  rate  or  loading  frequency,  specimen 
geometry,  and  flaw  size. 

Many  investigators  have  outlined  failure  models  for  advanced  fiber 
composites  considering  a  specific  composite  system,  a  specific  loading 
condition,  or  extension  of  metallic  (isotropic)  failure  theories. 
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An  outline  of  this  research  is  discussed. 


1 •  Theory  of  Failure 

Microscopic  fracture  analysis  has  been  the  subject  of  a  large  number 
of  publications  since  the  introduction  of  the  classical  theory  by  Griffith  (1). 
Irwin  (2),  Cherepanoe  (3),  and  Kostrov  and  Kikitin  (4)  improved  through 
their  analyses,  the  understanding  of  the  stress  state  in  the  vicinity  of  the 
crack.  Dynamic  crack  propagation  and  related  phenomena  have  been  reported 
by  such  researchers  as  Goldshtein  (5),  Cotterell  (6),  and  Daniel  (7). 

The  microscopic  stress  state  in  a  fibrous  composite  has  been  studied  by 
Koufopoulos  and  Theocaris  (8),  Leissa  et  al  (9),  and  MacLaughlin  (10). 

Research,  attempting  to  establish  bulk  composite  strength,  has  led  to  the 
publication  of  various  statistical  strength  models.  Two  of  these  are 
due  to  Scop  (11)  and  Friedman  (12).  Halpin  (13)  presented  an  empirical 
model  relating  failure  loads  and  time  effects  statistically. 

Macroscopic  failure  theories  have  been  developed  to  apply  to  a  given 
system  or  application.  In  general,  the  theories  of  failure  for  isotropic 
materials  have  been  extended  to  include  anisotropic  materials.  Most  research 
treats  the  composite  as  quasi-homogeneous  or  as  a  two-phase  system  (fiber 
and  matrix) .  Arbitrary  parameters  are  included  in  the  analyses  to  produce 
a  failure  surface.  The  Maximum  Distortional  Energy  Theory  (Von  Mises) 
is  the  basis  of  most  analyses.  Hill  (14)  extended  the  Von  Mises  theory 
to  include  anisotropic  materials.  Hill  produced  a  failure  surface.  Azzi 
and  Tsai  (15)  modified  the  Hill's  criterion  to  better  fit  test  data.  Chamis 
(16)  makes  several  assumptions  and  modifies  the  Von  Mises  theory  to  include 
both  tensile  and  compressive  loading  condition.  Tsai  and  Wu  (17)  considers 
tensile  and  compressive  loading  and  assumes  failure  will  exist  in  a  stress  space. 
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This  produces  a  modified  failure  surface.  These  models  are  essentially 
generalizations  of  the  Von  Mises  criterion  adapted  to  orthotropic  plates 
in  plane  stress.  The  models  do  not  specify  the  mode  of  failure  and  does 
not  distinguish  between  interfacial  bond  failure,  delamination  (fracture 
within  the  matrix),  and  simple  fiber  fracture. 

2.  Fracture  Mechanics  Applied  to  Composites 

Fracture  mechanics  was  developed  to  fill  the  need  within  isotropic 
pressure  vessels*  Wu  (18)  found  that  under  specific  conditions  isotropic 
fracture  mechanics  can  be  applied  to  anisotropic  materials.  These  conditions 
include:  (1)  definite  orientation  of  the  flaw  with  respect  to  the  principal 

asix,  (2)  stress  Intensity  factors  consistant  with  isotropic  materials, 
and  (3)  critical  crack  orientation  coincident  with  the  elastic  symmetry. 

These  conditions  satisfy  the  unidirectional  fiber  orientation  only. 
Unfortunately,  conventional  isotropci  fracture  mechanics  do  not  apply  in 
the  majority  of  anisotropic  structural  applications. 

Tetelman  (19)  outlined  the  potential  fracture  mechanisms  within  fiber 
composites  and  related  the  type  of  failure  to  its  constituent  characteristics. 
The  fiber,  matrix,  and  interfacial  behavior  under  load  generate  the  fracture 
condition.  Zweben  (2)  presents  an  analysis  of  composite  tensile  strength. 

Fibers  demonstrate  strength  and  length  effects  leading  to  at  least  three  modes 
of  composite  failure.  These  modes  are:  (1)  crack  propagation  in  the  matrix 
with  few  fiber  breaks,  (2)  fiber  break  propagation,  and  (3)  cumulative  damage. 

Waddoups,  Eisenmann,  and  Kaminski  (21)  apply  classical  fracture  mechanics 
to  laminated  composites.  A  circular  hole  with  flaws  transverse  to  the  principal 
axis  establishes  the  point  of  maximum  stress  concentration  and  defines  the 
crack  propagation  direction.  This  model  insures  fulfillment  of  the  specific 
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conditions  required  to  utilize  isotropic  fracture  mechanics.  The 
model  simplifies  the  analytical  approach  but  is  questionable  for  real 
application. 

Sih  et  al.  (22)  investigated  the  application  of  fracture  mechanics 
to  fiber  composite  systems  using  crack  models.  Existing  theory  agrees 
well  with  glass  fiber  composite  data  but  graphite  fiber  composites 
require  a  more  elaborate  model. 

Phillips  (23)  investigated  the  fracture  mechanics  of  carbon  fiber 
composites.  Crack  models  produced  stress  intensity  factors  which  agreed 
well  with  experimental  data.  Results  differed  using  natural  cracks  and 
"as-cut"  cracks. 

3.  Fatigue 

Fatigue  failure  in  metals  has  been  investigated  successfully  by  assuming 
that  a  portion  of  the  total  strain  energy  is  plastic  even  at  stresses 
below  the  proportional  limit.  Fatigue  failure  occurs  following  the  application 
of  cyclic  loads  over  a  period  of  time.  Saklind  (24)  reviewed  the  literature 
covering  fatigue  of  composites  and  existing  experimental  data.  Composites 
are  compared  to  metal  fatigue  criteria  and  composites  differ  in  low  cycle 
fatigue.  In  fact,  compsites  demonstrate  a  high  cycle  fatigue  strength 
greater  than  predicted  from  metals  criteria.  Low  cycle  fatigue  is  predicted 
as  a  critical  design  factor  for  composites. 

Holmes  and  Wright  (25)  tested  several  angle-ply  graphite- epoxy  specimens 
in  fatigue  and  concluded  that  the  observed  secondary  failure.  Primary 
failure  was  by  tension  at  the  reduced  cross-section. 

Several  strength  criteria  have  been  postulated  for  composites.  In  most 
cases  these  criteria  are  assumed  to  be  valid  for  static  as  well  as  dynamic 
loading.  Many  investigators  have  modified  the  Von  Mlses  yield  criterion  for 
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orthotropic  materials.  Sendeckyj  (26)  has  reviewed  these  and  other  strength 
criteria  for  composites. 

Hashin  and  Rotem  (27)  present  a  simple  failure  criterion  involving  a  fiber 
failure  mode  and  a  matrix  failure  mode.  The  failure  criterion  agrees  well 
with  glass  fiber  composite  fatigue  data. 
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SECTION  II 


TECHNICAL  APPROACH 

The  investigation  of  the  fatigue  strength  of  multilayer  advanced 
fiber  composites  should  consider  all  possible  variables.  These  variables 
include  composite  system,  fiber  orientation,  specimen  geometry,  hole 
size,  type  of  loading,  fatigue  loading  conditions  and  loading  frequency. 

To  observe  the  effect  of  individual  variables,  as  many  of  the  remaining 
parameters  must  be  held  constant  as  possible.  In  this  manner,  the  primary 
parameters  are  defined  and  incorporated  into  a  generalized  failure  theory. 

The  composite  system  considered  was  the  graphite^epoxy  (Thornel  300/ 
5208)  system.  Other  systems  should  be  investigated  to  establish  primary 
failure  modes. 

Two  test  specimen  configurations  were  used  in  the  research.  The 
static  specimen  was  6.5~inch  x  1-inch  and  eight  (8)  plies  in  thickness. 

These  specimens  contained  no  circular  holes.  The  fatigue  specimen  was  of 
similar  geometry  containing  a  3/16-inch  diameter  hole  in  the  center  position. 
Aluminum  grips  were  bonded  on  the  ends  of  the  specimens.  These  specimens 
are  shown  in  Figure  1. 

All  specimens  contained  a  general  fiber  orientation  of  [O/+0/O]^. 

The  angle  0  is  measured  between  the  longitudinal  axis  and  the  oriented 
fibers.  The  fiber  orientations  considered  were: 

[0]g  unidirectional 

[0/±15/0]^ 

[0/+30/0] 

—  s 

[0/+45/0]g 

[0/+60/0] 

—  s 


The  specimen  configuration  and  the  hole  size  are  fixed  with  varying 
fiber  orientations. 

The  static  test  specimens  (no  hole)  are  required  to  establish  the 
elastic  constants  necessary  for  analysis.  The  fatigue  specimens  (with  hole) 
are  used  to  investigate  the  fatigue  strength  and  failure  modes.  A 
sinusoidal  loading  function  of  basically  tension-zero'-tension  was  applied 
axially  to  the  specimens.  A  small  tensile  load  was  applied  at  the  minimum 
stress  to  insure  against  compressive  loads.  The  maximum  stresses  used  were 
70%,  62.5%,  and  55%  of  the  ultimate  tensile  stress  of  the  static  specimen. 
Two  loading  frequencies,  3  and  25  Hertz,  were  investigated  to  examine  low- 
frequency  and  high-frequency  tensile  fatigue  effects. 

The  static  and  fatigue  test  specimens  were  modeled  and  analyzed. 

Analytical  characterization  of  the  following  anisotropic  plates  were 

considered.  Included  were: 

static  -  no  hole 
dynamic  -  no  hole 

static  -  hole  infinite  plate 
dynamic  -  hole  infinite  plate 

static  -  hole  finite  plate 
dynamic  -  hole  finite  plate 

These  analytical  models  consider  all  possibilities  of  tensile  loading 
for  the  static  and  fatigue  test  specimens* 
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SECTION  III 


ANALYTICAL  CHARACTERIZATION 

The  analytical  characterization  of  the  load  and  deformation  distributions 

within  a  multilayer  fiber  composite  plate  was  examined  for  static  and 

dynamic  loading  conditions.  The  generalized  anisotropic  plate  analysis 

for  static  loading  conditions  without  a  flaw  (circular  hole)  was 

developed  to  define  the  field  equations  and  boundary  conditions  necessary 

for  further  analytical  development.  The  analysis  considers  the  linearly 

exact  relations  resulting  in  the  equilibrium  equations.  The  constitutive 

equations  are  outlined  for  a  laminated  medium.  The  field  equations  and 

boundary  conditions  are  based  on  the  Kirchhoff  hypothesis.  The  constitutive 

equation  coefficients  are  defined  for  the  composite  system  [O/+0/O]  . 

s 

This  composite  system  is  basic  to  this  research.  These  equations  are  basic 
for  the  static  loading  condition. 

The  field  equations  and  boundary  conditions  were  derived  for  a 
dynamically  loaded  anisotropic  plate  with  and  without  a  circular  hole.  The 
hole  in  an  infinite  plate  or  a  hole  in  a  finite  plate  becomes  a  matter  of 
definition  and  application  of  the  appropriate  boundary  conditions.  All 
systems,  whether  static  or  dynamic  loading  or  with  or  without  circular 
hole,  should  relate  to  the  generalized  anisotropic  plate  analysis  when  the 
proper  boundary  conditions  are  used. 

The  generalized  anisotropic  plate  analysis  defines  the  field  equations 
and  boundary  conditions  for  the  static  loaded  anisotropic  plate  (no  hole) . 

The  homogeneous  set  of  equations  was  analyzed  using  the  finite  difference 
method.  The  dynamic  stress  field  in  the  composite  considered  one-dimensional 
loading. 
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A  finite  difference  method  was  used  in  an  effort  to  solve  the  field 
equations  for  the  static  and  dynamic  models  containing  a  circular  hole. 
The  boundary  conditions  around  the  periphery  of  the  hole  were  difficult 
to  specify. 


^3  “  constant 
L  -  length 
M  -  stress  moment 


N  “  applied  moment 
P  ~  applied  tractions 
R  -  radius  of  hole 


t  -  time 


T  ~  stress  resultant 


u,v  ~  displacements  in  x  and  y  direction 
x,y,z  -  coordinates 
a, 3  -  constants 
c  ,  £  -  in-plane  strains 

XX  yy 

0)-  frequency 
p-  mass  density 
a-  applied  stress 
0 ,  (j)  -  angles 
V-  Poisson's  ratio 
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2.  Generalized  Anisotropic  Plate  Analysis 


2.1  Elastic  Constants 

The  generalized  anisotropic  plate  analysis  is  generated  for  a  plate 

of  thickness  **h"  and  the  reference  surface  is  at  the  mid  plane.  The 

linearly  exact  relations  are  defined  in  terms  of  stress  resultants  [28]. 

The  stress  resultants  are  defined  in  terms  of  stresses.  The  equilibrium  and 

compatibility  relations  are  defined  at  the  reference  surface.  Hooke* s  law 

is  written  for  an  orthotropic  material.  The  rotational  transformation  law 

for  elastic  constants  has  been  reported  by  Tsai  [29]. 

The  constitutive  equations  are  considered  for  a  layered  medium.  The 

field  equations  and  boundary  conditions  for  off  angle  layers  symmetrically 

layed“up  are  identified  from  the  analysis.  The  field  equations  can  be 

derived  for  the  static  as  well  as  the  dynamic  case. 

The  equations  for  the  determination  of  the  elastic  constants  by  axial 

test  can  now  be  considered.  For  purposes  of  determining  the  elastic  con- 
*  *  * 

stants,  b^^  ,  >  ^22  *  types  of  tensile  specimens  are  required. 

The  first  will  have  all  fibers  longitudinal  to  the  load  (0®)  and  the  other 
will  have  the  fibers  transverse  to  the  load  (90°).  A  sketch  of  the  coupon 
with  dimensions  and  coordinate  axis  is  shown  as  follows. 
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The  strain  expressions  which  satisfy  the  field  equations  and  boundary 


conditions  are  given  as  follows. 


22  £ 

(V22-S2'>  ' 


(2-1) 


.  =  12  ■  £ 

'yy  (c  c  -c  ^ 

^^11  22  12  ^ 


For  the  0°  direction 


*  :?c 

^^11  ^22  ~  ^12^ 


(2-2) 


:k  A 

^^11  ^22  "  ^12^ 


For  the  90®  direction 


(•’n  ”22  -  h2> 

-”12  ^ 

(”n”22-V  “ 


(2-3) 


Solving  for  elastic  constants,  the  result  becomes  the  following. 


e  °  e  90 
8bh  [  {-^)  ( — ^) 
^  P 


e  °  e  90  2 

4  90  ’  ^ 

P 


e  °  e  90 

8bh  [(-^)  (-2^9q) 

P 


^  e  “  e  90  2 

1  (_n.  +  _£Z  )  1 

4  ^  p°  90’  ^ 

P 


(2-4) 
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e  “  e  90 

_iz  +  yi— 

p°  90 
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e  °  e  90 
8bh  [(-^)  (-^2£_) 


1  E 

i  (jn. 

4  ^  p“ 


e  90 

+  _yz 

90 


)^] 


2.2  Governing  Equations  for  a  Laminate  with  Hole  and  Symmetrically  Loaded 

The  field  equations  and  boundary  conditions  are  outlined  for  a  dynamically 
loaded  8-ply  symmetrically  layed-up  composite  plate.  The  lay-up  for  the 
half  thickness  is  [O/+0/O].  A  circular  hole  of  radius  "R"  is  assumed  to 
exist  in  the  plate  and  a  uniaxial  membrane  excitation  is  imposed  on  one 
end  of  the  plate.  Thermal  effects,  bending,  body  forces  and  surface 
tractions  are  neglected.  A  diagram  of  the  plate  is  shown  as  follows. 


T 

L 

1 

J  j  1 1 

a(t) 


The  field  equations,  assuming  "equivalent  homogeniety"  through  the  thickness, 
become : 


,2 

c  ^  ^  f* 

^11  ^  2  ^66  ,,  2 

3x  by 


2  2  2 
9ul//--  I/-.  \  9v 

^^12  W"  "  ^ 


(2-5) 


^2  .2  .2  .2 
^^12  ^  ^66^  9x8y  ^  ^66g^2  ^22  ..  2  „.2 


9y 


9t 
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The  boundary  conditions  are: 
(1)  Stress  conditions  at  edges 


2,  t)  =  0 


Tj.^(R,  t)  =  0 


T  (x,+^,  t)  =  0 
xy  ’—2 


(R,  <}.,  t)  =  0 


(2-6) 


T  „  (±  y.  t)  =  0 


xy 


T  (±  y>  t)  =  8ha(t) 

XX 


(2)  The  stress  boundary  conditions  in  terms  of  displacements  are: 

'<=12  If  ^  <=22  If' 

(x,±  -5  t) 


+  f)  t  -  0 
’y  **  (K,±  t) 


3x  Bx"^  . ,  L  V 

(+  y,  t) 


(2-7) 


[C..  -^  +  ^ 

11  x  12  y  L 

(±  Y’  y»t) 


=  8h0(t) 


[(C^^  cos^<}.  +  C^2^in^^)  +  (2C^^cos<l>  sinc^  )  + 


o  o  3v 

(C22Sin  <p  +  37 


=  0 


(R,  <t>,t) 


[Cgg(cos^<J)  sin^<p)  1^  -  (C3^i-C^2^  ‘^°®***  ll 


C66(cos2^>  =  sin2<}>)  |j  +  (C22-C12)  cos(^sin(|r|^]  =  0 


Tho  bountlary  conditions  consider  a  free  outside  boundary  and  the  expressions 
are  the  same  wftli  or  without  a  circular  hole. 
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3»  Dynamic  Plate  Model  -  No  Hole 

The  homogeneous  set  of  equation  (2-5)  was  expressed  in  finite  difference 
form  and  was  solved  by  iteration.  The  finite  difference  form  of  equation 
(2-5)  is  as  follows: 

"3  "  ^u  ^1  u  -  2u  +  u 

^11  "  ^66^  ■ ;  f 


Ax 


Ay 


+  (C  +  C  )  (  ""6  ""lO  ~  "".12"  Is  ^ 
^  ^^12  ^66^  ^  4AxAy  ^ 


«=12  ^  S6>  <  -  ■  - - >  *  <=66  < - 

V  -  2v  +  V 

+  "22  '-^—2 — -  “ 

Ay 


%  ^1. 


Ax 


(3-1) 


where  the  general  fourth  order  finite  difference  representation  molecule  used 


for  internal  points  is  as  follows. 


Calculation  of  displacements  was  considered  for  the  upper-right  (+,+) 


quadrant  of  the  following  coupon  diagram. 


The  boundary  conditions  used  were: 

(1)  u  displacements  are  symmetric  across  x  axis  and  antis3nnmetric 
across  y  axis. 

(2)  V  displacements  are  S3mmietric  across  y  axis  and  antisymmetric  across 
X  axis. 

(3)  u  displacements  on  the  y  centerline  are  zero;  v  displacements  on  the 
X  centerline  and  at  x  =  L/2  are  zero. 

(4)  applied  displacement  at  x  =  L/2  is  u  =  u^  +  k  sinCojt 

where  is  initial  displacement,  =  0.01000  inches,  is  a 
suitable  constant  (0.00050  inches),  is  the  circular  forcing 
frequency,  t  is  time,  and  (j)-  is  a  phase  angle  (0  radians). 


(5)  the  free  boundary  of  the  quadrant  was  taken  as 
V  =  [k^  +  sin((0t  +  (j)2)]e“^  ~ 


(3-2) 


where  k^,  a»  3  are  constants  and  (j)^  is  a  phase  angle  (0  radians). 
This  boundary  condition  was  calculated  to  match  the  magnitude  of 
displacements  observed  at  the  y  centerline  of  the  static  tensile 
tests.  Displacement  u  on  the  free  boundary  was  determined  using 

euation  2-7.  The  equation  is: 

riu  iV]  ^ 

9x\Lb^. 

''2’2’ 


X  .1/2 


The  displacements  were  assumed  to  vary  as  from  a  value  of 

zero  at  X  “  0  to  a  maximum  given  by  the  boundary  conditions  evaluated 
at  X  =  L/2.  For  u  displacements,  the  following  equation  was  used. 


u  =  (k2  +  k^  sin(a)t  +  -e  )  + 

Uq  +  kj^  sinojt  + 


(3-3) 


Values  of  the  constants  used  were: 


-1.5000 


6.000 
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=  1.000 

=  0.5000 

(6)  the  V  displacements  along  the  y  centerline  were  taken  as  a 
linear  function  of  ”y"  ^.nd  the  u  displacements  were  assumed 
to  be  linear  functions  of  "x”  along  the  x  centerline.  Along 
the  X  centerline,  the  u  displacements  were 
U  =  ^  [Uq  +  sin(oi)t  +  (j)^^)] 

Along  the  y  centerline,  v  displacements  were: 

V  =  -  [k^  +  k^  sin(a)t  +  (1)2)]  e  ^ 


The  homogeneous  set  of  field  equations  corresponding  to  the  static 
case  was  solved  by  advancing  time  to  achieve  required  boundary  displacements 
and  then  iterating  for  the  displacement  field.  Calculations  were  made 
at  the  nodes  of  a  4x4  grid  in  the  upper  right  quadrant  (25  points). 
Convergence  giving  three  significant  figures  occurred  in  20  iterations. 
Approximately  10  seconds  of  computer  time  was  required  per  static  field. 

In  order  to  solve  the  nonhomo geneous  set  of  field  equations  (2-5) 
corresponding  to  time  dependency,  the  inertia  terms  of  the  right  hand 
side  of  the  equations  were  taken  as: 


8ph 


8ph 


2  =  f^(t)  f^(x,y) 

L 


2  2 
3  V  3  V 

By  solving  the  first  of  the  field  equations  for  and 


3  2 

X 


3  2 
y 


(3-4) 


and  sub¬ 


stituting  into  the  second  equation,  the  following  results. 


3  4 

X 


3^u 


^11  ^66 


2  2 
3x  3y 


<=11 
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(3-5) 


[  (- 

^11  ^66 


r  F,  (x,y))  + 

^12  ^  ^66  ^  ^  ^ 


2  2(F„(t)F  (x,y) 

<=,,  +  (Fl(t)F^(x,y)) - 1  -  0 

12  66 


'22 


3'u 


By  solving  the  second  equation  for  ^  and  ^and  substituting  into 

rv  ^  rv  ^ 


8x 


9y 


the  first,  the  following  results. 


SV  "66'  ("l2  "66>  -  "ll"22_ 

/.  I 


3  V 


^11^66 


2  2 
3x  3y 


(3-6) 


^22  9^  _  (^12/  ^66) 


C-,-,  a  ^ 

11  9y 


^11^66 


"l2  ^  "66 


'll 


^2  ''2<''’F)  -  ^  Fi(t)F,(x.y)l  -  0 


( 

^12  ■"  ^66  9y' 

The  following  substitutions  were  made  into  equations  3-5  and  3-6. 


Fi(t)  F^(x,y)  =  k^pho)'^  cos(a)^t)  sin(^  x)  cos(-^  y  -  x) 

n 

F„(t)  F-(x,y)  =  k  Phio^  cos(a)  t)  sin(^  cos(^  x) 


fiL  V  - 


(3-7) 


n 


n 


The  natural  frequency  of  the  coupon  is  and  is  equal  to  80- 00. 
Orthotropic  plate  velocity  was  used  by  Reuter  [30]  to  determine  displacement 
fields. 

Time  dependent  boundary  conditions  of  the  forcing  function  were  applied. 

A  time  increment  is  Ax  =  -r^  =  •  Addition  of  the  inertia 

2(0  ^ 

terms  required  further  iteration  to  establish  dispalcements  in  the  quadrant. 


Strain  energy  densities  were  calculated  for  the  one  quadrant. 
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The  strain  energy  density  fields  for  the  static  loading  (t=0)  of 
the  [0/^/0]  graphite-epoxy  composite  systems  are  shown  in  Figures  2-6. 

O 

The  graphite-epoxy  specimens  deform  approximately  10,000  y-inch  (at  x  = 

2  inch)  prior  to  fracture.  The  a  ,a  ,  and  a  stress  fields  were  also 
determined  but  not  included  in  this  report. 

The  dynamic  strain  energy  density  fields  for  the  graphite-epoxy 
composite  systems  were  determined  for  a  loading  frequency  of  50  Hertz, 
considering  time  steps  of  0.005,  0.010,  and  0.015  seconds.  The  static 


The  strain  energy  density  fields  for  the  three  (3)  time  steps  are  shown  in 
Figures  7-21.  The  change  in  the  strain  energy  density  contours  with  time 
notes  the  microscopic  variance  which  relates  to  fatigue  characteristics. 

A  failure  criterion  was  proposed  correlating  the  strain  energy  release 
rate  (G)  and  the  strain  energy  density.  It  is  assumed  that  the  strain  energy 
release  rate  is  proportional  to  the  strain  energy  density. 
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Graphite- Epoxy  Composite  under  10,000  y-inch 
Displacement  at  x  =  2-inch. 


0.500 


in-lb 


Epoxy  Composite  under  10,500  y-inch  Displacement 


0.500 


Figure  8.  Strain  Energy  Density  Field  of  Unidirectional 
Graphite-Epoxy  Composite  under  10,000  y-inch 
Displacement  at  x  =  2-inch 
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,500  y-inch  Displacement  at  x  =  2- inch 


Strain  Energy  Density  Field  of  [0/+15/0] 
Graphite-Epoxy  Composite  under  10,000 y -inch 
Displacement  at  x  =  2- inch. 
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Figure  12.  Strain  Energy  Density  Field  of  [0/±15/0]  Graphite-Epoxy 
Composite  under  9,500  y-inch  Displacement  at  x  =  2-lnch. 
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Figure  13.  Strain  Energy  Density  Field  of  [0/+30/0] 

Graphite-Epoxy  Composite  under  10,500 y -Inch 
Displacement  at  x  =  2-inch 
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Graphite-Epoxy  Composite  under  10,000y-inch 
Displacement  at  x  =  2-inch 


igure  15.  Strain  Energy  Density  Field  of  [0/±30/0]  Graphite-Epoxy 
Composite  under  9,500  p-inch  Displacement  at  x  =  2-inch. 


Composite  under  10,000  y-inch  Displacement  at  x  =  2-inch 


0.500 


Figure  18.  Strain  Energy  Density  Field  of  [0/145/0]^  Graphite-Epoxy 
Composite  under  9,500  y-inch  Displacement  at  x  •-  2~inch. 


er  10, 50' 


0,500 


Figure  20,  Strain  Energy  Density  Field  of  [0/±60/0]  Graphite-Epoxy 
Composite  under  10,000  y-inch  Displacement  at  x  ~  2-inch 


4.  Static  and  Dynamic  Plate  Model  With  Hole 

The  field  equations  (2-5)  were  converted  into  finite  difference  form. 
There  are  seven  (7)  groups  of  equations  that  must  be  considered  to  solve  the 
problem.  They  are: 

(1)  general  field  equations 

(2)  right  boundary  conditions 

(3)  top  boundary  conditions 

(4)  bottom  boundary  conditions 

(5)  left  boundary  conditions 

(6)  circular  hole  boundary  conditions 

(7)  points  have  close  relation  with  circular  boundary 

For  the  explicit  method,  is  required  to  be  small  enough  in  order  to 

generate  a  finite  difference  solution.  Time  steps  of  2x10  seconds  were 
required  to  achieve  a  reasonable  solution. 

The  grid  pattern  for  the  upper  right  quadrant  of  the  plate  is  shown 
in  Figure  22.  The  small  grid  around  the  hole  was  for  best  definition  of  the 
circular  boundary  conditions.  The  complication  of  defining  the  grid  points 
at  the  boundary  made  it  difficult  to  define  the  strain  energy  density  field 
or  stress  concentrations  around  the  hole. 

The  development  of  the  analysis  proved  frustrating  as  the  program 
generated  unreasonable  data  for  the  gird  pattern  and  boundary  conditions 
used.  Another  approach  is  being  investigated  to  minimize  the  problems 
generated  by  the  finite  difference  method. 
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Figure  22.  Grid  Pattern  for  Analysis  of  Composite  Plate  with  Hole. 


SECTION  IV 


I 


EXPERIMENTAL  DATA 

The  graphite-epoxy  system,  with  general  fiber  orientation  [O/+0/O]  , 

—  s 

was  subjected  to  tensile  static  fatigue.  The  test  specimen  geometries 

were  shown  in  Figure  1.  The  static  specimens  were  strain  gaged  with  EA-09- 

062TT-120  gages.  These  gages  provide  the  monitoring  of  both  longitudinal 

and  transverse  deformations.  The  strain  gaged  specimens  were  placed  in 

an  Instron  test  machine  and  deformed  at  0.02  inch/minute.  The  constant 

crosshead  movement  provides  determination  and  comparison  of  elastic  behavior. 

The  graphite-epoxy  fatigue  specimens,  with  3/16-inch  diameter  hole, 

were  subjected  to  tensile  cyclic  loading.  A  Gilmore  Universal  Testing 

System  was  used  to  apply  a  sinusoidal  load  function  to  the  specimens  under 

constant  load  conditions.  The  tensile  fatigue  tests  were  tension  (S  ) 

®  max 

-  0  -  tension  (S  ) .  The  minimum  tensile  stress  (S  .  )  was  not  zero  but 
max  ^  min 

slightly  in  tension  to  avoid  microbuckling  from  compressive  loading. 

The  value  of  R  was  equal  to  0.005  in  all  tests  (R  =  S  .  /S  ).  The  fatigue 

min  max 

specimens  were  subjected  to  3  loading  conditions  (55%, 62. 5%,  and  70%  of 
static  ultimate  tensile  stress)  and  2  frequencies  (3  and  25  Hertz).  The 
loading  conditions  were  based  on  static  specimen  geometry.  The  3/16-inch 
diameter  holes  were  not  considered  when  calculating  the  cross-sectional  areas. 
The  specimens  were  allowed  to  be  tested  until  failure  occured  or  10  cycles 
were  attained.  The  tensile  fatigue  data  provided  an  insight  to  the  failure 
mechanism. 

The  test  specimens,  tensile  static  and  fatigue,  were  examined  following 
fracture  in  order  to  note  the  failure  surfaces. 
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1,  Graphite-Epoxy  Composite 


The  tensile  static  and  dynamic  test  specimens  were  supplied  by  the 
Air  Force  Flight  Dynamics  Laboratory.  Thornel  300  graphite  yarn  was  the 
fiber  combined  with  5208  epoxy  resin.  The  properties  of  Thornell  30  [31] 
are  contained  in  Table  1.  The  Thornel  300/5208  composite  flat  plate 

Table  1.  Properties  of  Thornel  300  Graphite  Yarn 


Property 

Unit 

Value 

tensile  strength 

lb/ in^ 

360,000 

tensile  modulus 

Ib/in^ 

34  X  10^ 

density 

g/cc 

1.76 

elongation  at  break 

% 

1 

specimens  were  autoclave  cured  (Whittacker  Corp.  specification)  and  post 
cured  for  4  hours  at  400°F.  The  individual  constituents  within  the  laminates 
are  summarized  in  Table  2. 

Table  2.  Content  of  Thornel  300/5208  Laminates 


Fiber 

Orientation 

Specimen 

Letter 

v/O  Fiber 

v/O  Resin 

v/O  Void 

[0/+0/0]^ 

A 

68.3 

30.4 

1.3 

[0/+15/0] 

—  S 

B 

65.1 

33.8 

1.1 

[0/+30/0]g 

G 

66.8 

32.1 

1.1 

[0/+45/0]g 

D 

66.9 

31.5 

1.6 

[0/+60/0] 

““  S 

E 

65.2 

33.3 

1.5 

The  laminates  vary  from  65.1  to  68.3  in  fiber  volume  percent. 

This  variation  is  small  so  the  test  results  of  the  laminates  can  be 
compared . 

The  fabrication  of  theThornel  300/5208  laminates  was  of  high  quality; 
which  is  basic  to  a  test  program. 
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2. 


Static  Test  Results 


Two  each  of  the  [0]q  (A)  and  [0/+15/0]  (B)  static  test  specimens 

o  —  s 

were  subjected  to  uniaxial  tensile  loading.  Strain  gages  were  placed  on 
each  flat  surface  of  the  test  specimen.  Differentials  in  strain  gage 
readings  would  indicated  bending  in  the  specimen.  The  average  of  the  two 
gages  on  each  specimen  was  used  to  produce  the  stress-strain  curve. 

The  stress-strain  relationships  for  these  fiber  orientations  are  shown  in 
Figures  23  and  24. 

Two  each  of  the  [0/+30/0]  (C) ,  [0/+45/0]  (D) ,  and  [0/+60/0] 

s  s  s 

(E)  static  test  specimens  were  subjected  to  uniaxial  tensile  loading. 

One  strain  gage  was  placed  on  the  flat  surface  of  each  test  specimen  for 
monitoring  longitudinal  and  transverse  deformations.  The  stress-strain 
relationships  for  these  fiber  orientations  are  shown  in  Figures  25,  26  and 
27. 

The  mechanical  properties  of  the  [0/+9/0]^  Thornel  300/5208  composite 
specimens  are  summarized  in  Table  3.  The  data  recorded  are 


Table  3.  Mechanical  Properties  of  Thornel  300/5208  Composite  Specimens 


Property 

[0]„ 

A  ® 

[0/+15/0]^ 
B  ® 

[0/+30/0] 

C  ® 

[0/+45/0]^ 

D  ® 

[0/+60/0] 

E  ® 

^^Jvs±) 

227,000 

155,480 

139,910 

116,790 

100,360 

£^^(y-in/in) 

9830 

7711 

9610 

8710 

8540 

E^(10^  psi) 

21.9 

19.3 

15.4 

12.5 

11.8 

V 

0.29 

0.64 

0.98 

0.65 

0.27 

Fpj^(psi) 

96,500 

82,000 

the  average  of 

two  specimens 

per  fiber 

orientation. 

The  mechanical  properties  generated  with  the  static  test  specimens 
were  used  in  tensile  fatigue  testing. 

3.  Tensile  Fatigue  Results 


Tensile  fatigue  data  were  generated  for  the  Thornel  300/5208  composites 

using  a  sinusoidal  loading  function  under  constant  load  conditions. 

Three  loading  conditions  and  two  frequencies  were  used.  The  S  used 

max 

was  55%,  62.5%,  and  70%  of  tensile  ultimate  stress  determined  by  static 
test  and  these  static  test  specimens  were  used  as  control  specimens. 
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Stress  (ksi) 


Figure  24.  Stress-Strain  Curve  for  [Olg  Thornel 

300/5208  Composite 


46 


Stress  (Ksi) 


Figure  24.  Stress-Strain  Curve  for  [0/+15/ 
Comoosite. 
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Longitudinal  Strain  (y-in/in  x  10  ) 

re  29.  Stress-Strain  Curve  for  [07+30/0]^  Thornel  300/5208 
Composite. 
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Stress  (ksi) 


Figure  26..  Stress-Strain  Curve  for  [0/+45/0]  Thornel  300/5208 
Composite. 
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The  applied  stress  considered  the  plate  withoul  a  circular  hole.  The 
maximutn  alternating  stresses  for  the  various  geometries  are  given  in 
Table  4.  The  stresses  were  obtained  from  no  hole  specimens  and  applied 
to  specimens  with  hole.  The  value  of  R  was  equal  to  0.005.  The  specimens 
were  submitted  to 

Table  4.  Maximum  Alternating  Stresses  Used  in  Tensile  Fatigue  Testing  (No  Hole) 


Property 

[0/+0/0]^ 

A  ® 

[0/+15/0]^ 
B  ® 

[0/+30/0]^ 

C  ® 

[0/+45/0] 

D  ® 

[0/+60/0]^ 

E  ® 

F^^(psl) 

227,000 

155,480 

139,910 

116,790 

100,360 

70%F^ 

tu 

157,000 

108,000 

96,500 

82,000 

68,000 

62.5%F^ 

tu 

141,300 

97,200 

86,850 

73,800 

61,200 

55%F 

tu 

125,600 

86,400 

77,200 

65,600 

54,400 

constant 

load  conditions  at  3  and 

25  Hertz.  The 

results  of 

the  tensile 

fatigue  testing  are  shown  in  Figures  28  through  32.  The  ultimate 

tensile  strength  of  the  static  specimen  without  holes  is  noted  on  the 

6 

plots.  The  test  was  completed  when  either  failure  occurred  or  10  cycles 
were  obtained.  Two  (2)  specimens  were  tested  per  condition.  As  a  result, 
no  statistical  evaluation  of  the  data  is  possible.  However,  the  limited 
data  do  indicate  trends. 

The  damage  indicated  by  cyclic  loading,  as  compared  to  static  loading 
is  less  for  the  D  and  E  configurations.  This  relates  directly  to  the 
stress  concentrations  present  at  the  periphery  of  the  hole.  The  life  (number 
of  cycles  to  failure)  was  less  for  low-frequency  fatigue  specimens  (3  Hertz) 
than  the  high-frequency  fatigue  specimens  (25  Hertz).  All  high-frequency 
fatigue  specimens  attained  10^  cycles  with  a  of  low-frequency 
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[0/+15/0].  Thornel  300/5208 


(Ls>|)  ss0aq.s  aii-suai  uinuiLXBW 
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Figure  30.  Tensile  Fatigue  Curve  for  [0/+30/0]  Thorne!  300/5208  Composite. 


[0/+45/0]^  Thorne!  300/5208 


fatigue  specimens  all  failed,  even  under  the  low  S 

max 

The  damage  mechanisms  occurring  within  the  composites  appear  to  be 
related  to  frequency  (time)  and  level  of  loading.  However,  the  mechanisms 
appear  to  be  of  similar  type, 

4,  Test  Specimen  Failure 

The  tensile  static  and  fatigue  test  specimens  were  observed  during 
test  and  following  fracture  to  note  mode  of  failure.  Test  specimen 
notation  was  used  for  identification  by  AFFDL.  Fractured  static  test 
specimens  are  shown  in  Figures  33  through  35.  The  tensile  fatigue 
specimens  are  shown  in  Figures  36-45. 

The  static  specimens  reached  a  maximum  load  or  deformation  and 
fractured  suddenly.  The  unidirectional  specimens  splintered  while  fibers 
broke.  The  B  specimens  fractured  along  the  +15®  fiber  orientation. 

The  specimens  containing  +30®  and  +45®  plies  fractured  almost  in  a  straight 
line  across  the  width.  The  specimen  with  +60®  plies  fractured  near  the  grip  at 
the  +60®  angle. 

The  tensile  fatigue  specimens  exhibited  secondary  as  well  as 
primary  modes  of  failure.  Fiber  breakage  was  observed  in  all  tensile 
fatigue  tests.  Separation  of  portions  of  the  specimen  from  the  grips 
was  noted  in  many  cases.  Delamination  of  portions  of  the  surface 
plies  was  observed  in  C,D,  and  E  specimens.  This  is  related  to  Poisson’s 
effect.  Fiber  breakage  and  delamination  was  observed  on  the  edge  of 
all  specimens.  All  of  these  failure  modes  are  considered  secondary. 
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Figure  35.  Static  Test  Specimens: 

(a)  Fatigue  Specimen,  (b)  A-1-26  (with  hoi 
(c)  B-1-2  (with  hole). 
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Figure  36.  Tensile  Fatigue  Specimens  A  (25  Hertz); 

(a)  70%  F^^^,  (b)  62.5%  F^^,  (c)  55%  F^^. 


Figure  37, 


Tensile  Fatigue  Specimens  A  (3  Hertz); 
(a)  70%  F*. (b)  62.5%  F^.  .  (c)  55%  F^ 


Figure  38.  Tensile  Fatigue  Specimens  B  (25  Hertz); 

(a)  70%  F^jj,  (b)  62.5%  F^^^.  (c)  55%  F^^^ 


Figure  39,  Tensile  Fatigue  Specimens  B  (3  Hertz); 

(a)  70%  F.  ,  (b)  62.5%  F.„,  (c)  55%  F^ 


Figure  42.  Tensile  Fatigue  Specimens  D  (25  Hertz) 
(a)  70%  F.  ,,  (b)  62.5%  F.  .  (c)  55%  F. 


Figure  43.  Tensile  Fatigue  Specimens  D  (3  Hertz) 
(a)  70%  F.  .  (b)  62.5%  F., ,,  (c)  55%  F 


SECTION  V 


DISCUSSION 

The  analytical  characterization  of  the  load  and  deformation 
distributions  within  a  multilayer  fiber  composite  plate,  with  and  without 
a  circular  hole,  under  both  static  and  dynamic  loading  conditions  was 
initiated.  The  generalized  anisotropic  plate  analysis  for  static 
loading  conditions  without  flaw  (circular  hole)  was  developed  to  derive 
the  field  equations  and  boundary  conditions  necessary  forfurther 
analytical  development.  This  generalized  anisotropic  plate  analysis 
has  been  reported  in  various  forms  in  the  technical  literature. 

The  field  equations  derived  were  used  in  the  analytical  development 
of  the  dynamic  plate  model  (no  hole).  The  field  equations  were  expressed 
in  finite  difference  form.  A  homogeneous  set  of  field  equations  corresponding 
to  the  static  case  was  solved  by  incrementing  time  to  achieve  required 
boundary  displacements  and  then  iterating  for  the  displacement  field-  The 
strain  energy  density  fields  for  both  static  and  dynamic  loadings  of 
a  unidirectional  graphite-epoxy  were  developed  for  several  initial  displacements. 
The  correlation  of  strain  energy  density  with  the  strain  energy  release  rate 
(G)  would  be  a  basis  for  a  theory  of  failure  in  a  future  study. 

The  generalized  anisotropic  plate  analysis  produced  field  equations 
which  were  converted  into  finite  difference  form  in  order  to  develop 
the  static  and  dynamic  plate  model  with  a  circular  hole.  Specification 
of  the  boundary  conditions  at  the  circular  hole  was  not  properly  modeled. 

As  a  result,  the  strain  energy  density  fields  and  stress  concentrations 
around  the  hole  were  not  calculated. 

The  Thornel  300/5208  static  and  tensile  fatigue  test  specimens  were 
provided  by  the  Air  Force  Flight  Dynamics  Laboratory.  Sendeckyj ,  et  al 


65 


[32]  stated  that  the  tensile  modulus  of  the  unidirectional  composite 
was  23-25x10^  psi.  This  was  slightly  higher  than  21.9x10  psi.  Also 
noted  were  a  shear  modulus  (G)  of  0.77x10^  psi  and  tensile  modulus  of 
90“  fiber  oriented  composite  of  1.6x10^  psi.  The  fiber  orientations 
containing  +30“  and  +45“  plies  demonstrated  tensile  proportional  limits. 

The  other  static  specimens  demonstrate  a  stiffening  effect  rather  than 
a  proportional  limit. 

The  tensile  fatigue  data  were  generated  for  five  fiber  orientations, 

three  loading  conditions,  and  two  frequencies.  The  was  55%, 

62.5%,  and  70%  of  tensile  ultimate  stress  determined  by  static  test. 

The  applied  stress  considered  the  plate  without  a  circular  hole.  For 

all  fiber  orientations  tested  at  25  Hertz,  the  55%  specimen  did  not 

fail  at  10^  cycles.  For  70%  F^.^,  the  specimens  failed  after  a  few 

thousand  cycles.  Testing  at  3  Hertz  produced  lower  fatigue  lines.  At 

70%  F  ,  failure  occurred  after  a  few  hundred  cycles  for  all  fiber  orientations 
tu 

The  55%  F  specimens  failed  between  200,000  to  300,000  cycles  except 
tu 

[0/+45/0]g  which  failed  around  450,000  cycles.  The  low  frequency  fatigue 
specimens  have  lower  lives  than  the  high-frequency  fatigue  specimens. 

The  damage  mechanisms  occurring  within  the  composites  appear  to  be  related 
to  frequency  (time)  and  level  of  loading.  The  mechanisms  appear  to  be  of 
a  similar  type. 

The  unidirectional  (A)  Thornel  300/5208  composites  failed  under 
maximum  strain  conditions  in  static  and  tensile  fatigue  loadings. 

Secondary  effects  included  fiber  breakage  and  pullout  from  the  grips. 

The  degree  of  fragmentation  was  related  to  loading  conditions  and  frequency. 
The  static  specimen  the  least  fragmentation,  followed  by  low-cycle  then 
high-cycle  fatigue  specimens.  Static  specimen  C  failed  by  the  +30“  plies 
fracturing  causing  the  0“  plies  to  fail.  The  tensile  fatigue  specimens 
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failed  in  a  similar  manner  with  the  crack  originating  at  the  center  hole. 

Static  specimen  D  failed  cleanly  across  the  width.  Secondary  failure  effects 
were  noted  in  the  tensile  fatigue  specimens.  These  included  delamination 
(due  to  Poissons  effect),  fiber  breakage,  and  edge  fiber  delamination. 

The  crack  was  initiated  at  the  center  hole  and  propagated  between  the  +45 
and  the  -45®  plies  longitudinally  until  fracture  occurred.  Static  specimen  E 
failed  in  the  +60®  direction.  The  tensile  fatigue  specimens  failed  similar 
to  the  D  specimens  with  the  same  secondary  effects  and  longitudinal  delamination. 

This  research  has  outlined  several  variables  which  influence  the 
failure  of  graphite-epoxy  composites.  The  data  indicates  trends  not 
specifics  to  establish  a  definite  theory  of  failure. 
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SECTION  VI 


CONCLUSIONS 

The  analytical  characterization  of  the  load  and  deformation  distributions 
within  a  multilayer  fiber  composite  plate  (without  hole)  under  both  static 
and  dynamic  loading  conditions  was  accomplished  using  the  finite  difference 
method. 

Static  and  tensile  fatigue  testing  of  Thornel  300/5208  composites 
with  and  without  holes,  was  accomplished.  Tensile  fatigue  testing  involved 
five  fiber  orientations,  three  loading  conditions,  and  two  frequencies.  Speci¬ 
mens  tested  at  3  Hertz  had  lower  tensile  fatigue  properties  than  those 
tested  at  25  Hertz. 

The  failure  surfaces  of  the  specimens  demonstrated  the  effect  of 
testing  conditions.  Secondary  failure  mechanisms  such  as:  delamination, 
fiber  breakage,  and  edge  fiber  delamination,  were  present  in  D  and  E  specimens. 
Longitudinal  delamination  of  the  D  and  E  specimens  was  unexpected.  The 
static,  low-frequency  fatigue,  and  high  frequency  fatigue  failures  were 
generally  similar,  differing  only  in  degree. 
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